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CN| , The ability to now make measurements of Be and B as well as put 

^ I constraints on ^Li abundances in metal-poor stars has led to a detailed 

1-^ ■ reexamination of Big Bang Nucleosynthesis in the A^ 6 regime. The 

nuclear reaction network has been significantly expanded with many 
new rates added. It is demonstrated that although a number of A > 7 



c/3 ! reaction rates are poorly determined, even with extreme values chosen, 

^ ' the standard homogeneous model is unable to produce significant yields 

(Be/H and B/H < 10~^^ when A < 7 abundances fit) above A — 7 
^ ■ and the ^Li/^Li ratio always exceeds 500. We also preliminarily explore 

inhomogeneous models, such as those inspired by a first order quark- 
hadron phase transition, where regions with high neutron/proton ratios 
can allow some leakage up to A > 7. However models that fit the A < 7 
abundances still seem to have difficulty in obtaining significant A > 7 
yields. 



1. Introduction 

Over the last quarter century the standard homogeneous model of Big Bang nucleosyn- 
thesis has proved spectacularly successful at predicting the primordial abundances of the 
light elements. In particular, Homogeneous Big Bang Nucleosynthesis (BBN) (see Walker 
et al. 1991 and references therein) successfully fits ^H, ^H, ^H, ^He, and ^Li abundances 
in primordial objects (extremely low Z) over a dynamic range in abundance of almost ten 
orders of magnitude. The fit to these abundances has become the prime determinator of 
the cosmological density of baryons (O^ ~ 0.05 where Qb is the fraction of the critical den- 
sity in baryons). Also of significance has been BBN's prediction of the number of neutrino 
fiavors (Steigman, Schramm & Gunn 1977, Yang et al. 1984, Walker et al. 1991) which has 
now been confirmed in accelerator experiments. 

Because of the double particle instability gaps at A = 5 and A = 8 it was rapidly 
recognized that homogeneous BBN produces small amounts of A = 7 (^Li/H~ lO"^'^) and 
no significant yields for A > 7. Recently observers have begun to be able to observe Be 
and B in extreme pop. II stars with low Z (Rebolo et al. 1988, Ryan et al. 1990, Gilmore, 
Edvardsson & Nissen 1991, Duncan, Lambert & Lemke 1992, Gilmore et al. 1992, Ryan et 
al. 1992) and have begun to put limits on ^Li in similar objects (Andersen, Gustafsson & 
Lambert 1984, Spite & Spite 1982, Hobbs & Pilachowski 1988, Hobbs & Thornburn 1991). 
Data to date on Be and B are best understood as being due to cosmic ray spallation 
(Walker et al. 1992, Steigman & Walker 1992) and ^Li is still undetected. However, 
because of the potential for new developments here, we have reexamined BBN yields with 
a particular focus on A > 6. While the basic conclusions of the earlier calculations remain 
unchanged, we did note that the networks used earlier were not particularly complete for 
A> 6. (This relative incompleteness had little or no effect on A < 6 yields.) Therefore we 
have extended the nuclear reaction network and included many links that were missing (or 

2 



poorly estimated) in earlier calculations. Where links in our new network are poorly known 
we have run the calculation with a range of values to assess the sensitvity to those links. 
We feel this present paper is the most thorough exploration of high A{> 6) Homogeneous 
Big Bang Nucleosynthesis done to date. 

A much discussed alternative to Homogeneous Big Bang Nucleosynthesis has been the 
first-order quark-hadron phase-transition inspired inhomogeneous model (Alcock, Fuller & 
Mathews 1987, Applegate, Hogan & Scherrer 1988, Turner 1988, Terasawa & Sato 1989, 
Kurki-Suonio et al. 1990). It had been proposed (Boyd & Kajino 1989, Malaney & Fowler 
1989) that in such a model, the high n/p regions may allow leakage beyond the A = 8 
gap and produce interesting amounts of Be and B. While Terasawa & Sato (1990) argued 
against such leakage, we noted that their network was not as complete as the one we 
have developed for the homogeneous case and that our more complete network may be 
needed to fully explore the situation. In addition, with our enhanced network, we are 
able to quantitatively compare the yields of the heavier elements in the homogeneous and 
inhomogeneous models. 

Thus in this paper we will also do a preliminary exploration of high n/p conditions 
with a homogeneous calculation. (Since high n/p values only occur in the low density 
zones of inhomogeneous models, high baryon to photon results are less significant than the 
low baryon to photon results for our high n/p calculation.) In a future paper we will do a 
more complete exploration of inhomogeneous models using a full multizoned model with 
back diffusion as in Alcock, Fuller & Mathews (1987), Kurki-Suonio et al. (1990), Terasawa 
& Sato (1990). However, at present to get a feel for the full network and see the critical 
points in the calculation we felt this simpler exploration was necessary and illustrative of 
the network itself. 

One point we emphasize is that the high A calculations cannot be used without making 



sure that the low A abundances are fit. As Kurki-Suonio et al noted, the inhomogeneous 
calculations still require about the same O^ as the standard homogeneous BBN model. 
Thus any leakage up to high A must occur with standard Qb values if it is to be relevant. 

2. The Reaction Network 

Up until now, the emphasis in measuring astrophysically relevant nuclear reaction 
rates has generally been on those reactions involved in the production of elements up to 
Li. Of the yields for A < 7 only ^Li has been found to be sensitive (Kawano et al. 1988, 
Krauss & Romanelli. 1990) to uncertainties in the measured reaction rates. 

We have updated the reaction network in an attempt to resolve this question for 
A > 6. In table 1 we list the reactions that have been added (A) or updated (U) to 
the last version of the code (Walker et al. 1991) which itself is an updated version of 
Wagoner's code (Wagoner 1969), modified to evaluate neutron-proton weak interaction 
rates according to the prescription of Walker et al. (1991). The majority of the reaction 
rates used came from Caughlan and Fowler (1988). To this set we added reactions from 
Wagoner (1969), Lederer & Shirley (1978), Endt & Van der Leun (1978), Ajzenberg-Selove 
(1983), Tuh (1985), Malaney & Fowler (1989), Boyd & Kajino (1989), Wiescher, Steininger 
& Kappeler (1989), Wang, Vogelaar & Kavanagh (1991), Kawano et al. (1991), Brune, 
Kavanagh, Kellogg & Wang (1991), Barhoumi et al. (1991), Becchetti et al. (1992), Boyd 
et al. (1992a), and Boyd et al. (1992b). We also included reactions from Smith, Kawano & 
Malaney (1992) as given in Kawano (1992). In a number of cases (listed in table 1) we were 
unable to find any reliable figures for reactions in the relevant temperature range and made 
estimates (E). For a few reactions we have found more than one rate published recently 
and have tested for sensitivity (S) when these differed beyond stated errors. The complete 
network is shown in fig. 1. For the most part we used the most recently published rate, 
whenever more than one expression was available. For the reactions labelled S in table 1 



we ran the code using the rates from each of the sources hsted. For estimated rates or 
rates where experimental values differed outside of the published error bars we made tests 
of the sensitivity of our results by arbitrarily increasing and decreasing that particular rate 
by a factor of 1000 or by using the extreme experimental value to see the sensitivity of the 
result. The actual prefered reaction rates selected are shown in table 1, and the explicit 
rates used are given in Appendix 1. 

In order to ensure that our network is sufficiently extensive, we have plotted flow 
diagrams (see figs. 2) for rjio = 1.0 for both the standard calculation, and for high n/p 
(^10 = ?7 X lO^'^). These were produced by connecting the nuclide with the greatest increase 
in mass fraction, to that with the greatest decrease (omitting nuclides with A < 4) at each 
time step. This gives a pictorial representation of which links are most significant in 
producing any given nuclide. Fig. 2a shows that in the standard model, most of the fiow 
proceeds along the central portion of the network, with the exception of the fiow to ^^O. It 
was this 1^0 that prompted us to add ^^C, ^^N, ^^N, ^^O, ^OQ, ^^F, ^Op, sip and ^iNe to 
our original network, however this made no difference to the yields of the light elements. 
Fig. 2b shows that for n/p = 10, even with the extended network, there is a significant 
fiow on the neutron-rich side, particularly to ^^O and ^^O. Obviously, to accurately explore 
neutron-rich fiows requires networks as rich as the one we use here. Previous explorations 
of A > 7 yields have not used such an extensive network. 

3. The Results 
The yields of the light elements are shown as functions of the baryon to photon ratio 
?7, for a neutron mean-life of 889.6 sec. in figs. 3. Figure 3a shows the "^He mass fraction as 
a function of the baryon to photon ratio, t] = ub/ti-^ and figure 3b shows the abundances 
of ^H, ^He, ^Li and ^Li. In particular, ^Li is always greater than ^Li by a factor of at 
least 500 (at 7710 = 0.01), while for 7710 ~ 3 (where calculations agree with observations) 
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the ^Li/^Li ratio is 4700. 

Of the reactions labelled S in table 1, the only rate whose variation has a significant 
effect on the results is ^Li(t,n)^Be. This sensitivity is shown in fig. 3c. Furthermore, 
variation in this reaction affects only ^Be and ^°B. Yields of the lighter elements (H, He, 
Li) are unaffected by the variation in any of the rates S. 

The double hump for ^^B is a result of the fact that it can be created directly (for low 
1]) or through ^^C (high i]) which then /3-decays to ^^B. The highest yield for both ^Be and 
^°B is given by the ^Li(t,n)^Be rate taken from Boyd & Kajino (1989). The lowest is from 
Malaney & Fowler (1989). Note that in no case is it possible to produce a ^Be number 
density relative to hydrogen greater than about 10~^^. If we consider the observational 
limits on H, He and Li then ry is constrained by 2.8^ ^lo^ 4.0 (Walker et al. 1991). In 
this range ^Be/H has a maximum yield of 6 x 10~^^, and ^^B/H a maximum of 8 x 10~^^. 

For the reactions labelled E it was necessary to make estimates of the rates. Unfor- 
tunately the presence of resonances at nucleosynthesis energies can influence the reaction 
rate by many orders of magnitude. We have tried a number of approaches to this problem. 
Some of the reactions in question are given in Delano (1969). However, Delano parame- 
terizes his rates for 1 < Tg < 10 and is not explicit about his estimation methods. Nor are 
they transparent from his rate expressions which are all fit to a standard form. The only 
theory he mentions is a statistical nuclear model which is only accurate for heavy nuclei 
with a large, even distribution of levels lying within nucleosynthesis energies. 

To understand the possible effects of the unmeasured rates, one would ideally like to 
at least have upper bounds on them. Such limits are enormous, however, as resonances at 
the effective nucleosynthesis energy 

E,s=[naZoZ,kT{ficy2y/'f/' = 0.12204Z^^/'z^,/ \'^'^\ T^/^MeV (1) 

(where a ~ 1/137, Zq and Zi are the charges of the incoming nuclei, and fx their reduced 
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mass) can enhance a reaction rate to many orders of magnitude larger than reactions 
that are nonresonant, or occur through the tail of a resonance. We have made estimates 
of the rates to compare with Delano's. To do this we assumed each reaction had two 
contributions, one nonresonant and one resonant. In both cases the reaction rates can be 
given by expressions of the form 

(at;)T = AoT-'^exp(-a/T"^) (2) 

where the constants Aq, a, n and m differ for the two cases. 

In the nonresonant case, the exponential part of the rate is determined by the nuclear 
masses and charges, so that it is only the preexponential factor which is unknown. We 
have estimated this factor (the unknown part of which is the cross section factor 5'(0)) 
using the method of Fowler & Hoyle (1964), and making conservative estimates of optical 
model parameters we have 

(^av) ~ - 777; 7T7T7TT79 exp(2a: - r) (3) 



,//; \kTj {Woy/^ 

where Vq ~ 40 MeV gives a good fit to known cross-sections 



x = 2\/2fiZoZie^R/n'^; R = Ro{AI^^ + aY^); Rq = lAMm (4) 



and 

r = 3E,ff/kT (5) 

Thus our rate is computed in an appropriate low-energy limit, but arises from an optical 
model which assumes a smoothly varying and large level density which we do not expect 
to be very accurate for the light nuclei we consider. 



In the resonant case, we used a recent tabulation by Ajzenberg-Selove (1985, 1986, 
1988) to locate levels and to find their widths. Knowing the resonant energy fixes the 
exponential term in the rate, and so again we are uncertain only in the preexponential 
factor. The unknown part here is the reduced width (0^7) r ~ arErVr, in which the cross 
section at resonance, a^ is unknown. For this we have assumed an arbitrary but generous 
value of 1 barn. 

Using these estimates can have a significant effect on nuclear yields (although not for 
^Be or any of the lighter elements). In particular, ^'^B(q;,7)^'^N and ^^B(a,7)^^N reduce 
the yields of ^*^B and ^^B respectively by two or three orders of magnitude at 7710 = 3, 
and while ^^C(q;,7)^^0 has no effect at ryio = 3, it removes ^^C almost completely above 
ryio = 10. We emphasize however, that these reaction rates are dominated by contributions 
from resonances near the entrance channel, and that we consider them to be extreme upper 
limits. 

In order to find more reasonable limits we have also tried using Delano's expressions 
directly (with results identical to those obtained by omitting all reactions E) and estimating 
rates on the basis of "similar" reactions. For each reaction E, we used a rate which was 
higher than the rate of any reaction with a similar form (say all reactions X(n,p)Y for the 
case of ^^N(n,p)^^C), and also higher than Delano's rate (where one was available) within 
the range 1 < Tg < 10, and than increased the result by a factor of 1000. The only reaction 
for which this had any effect on the results was ^Li(d,n)^'^Be. Using 1000 times the rate 
for ^Li(d,n)^Be (Malaney & Fowler 1989) increases the yield of ^*^B by up to a factor of 2 
for ?7io^ 0.1. (For ?7io^ 0.1 the effect is again insignificant.) 

Clearly further experimental data would be helpful, however we believe (with a few 
special exceptions that will be mentioned later) that these reactions are likely to have little, 
if any, effect on the light element abundances. 



The significance of tfie Be and B calculations has been made evident by recent obser- 
vations in a number of low metallicity halo dwarf stars of these two elements. It is generally 
known that big bang nucleosynthesis is incapable of producing an observable amount of 
either of these two light elements. In fact, Be and B have generally been thought to have 
been produced by cosmic ray spallation. Indeed, in recent analyses, the observed Be and 
B have been argued to be explicable entirely in terms of cosmic ray spallation in the early 
galaxy, while maintaining consistency with big bang nucleosynthesis . (Spallation also pro- 
duces ^Li which thereby reduces the required production from big bang nucleosynthesis, 
but as emphasized by Olive & Schramm (1992) this reduction is still completely compatible 
with the other cosmological light element abundances). Furthermore, it has been specu- 
lated that inhomogeneous models may provide for enhanced Be and B abundances relative 
to the standard model. Here, we have found that the standard model production of Be 
and B is indeed negligible relative to the observations and in our exploration of high n/p 
we find that while Be and B yields are enhanced, they are still well below the observations. 

4. Limits on Inhomogeneous Yields 

To obtain an extreme upper limit on the yields produced by inhomogeneous nucleosyn- 
thesis, we show in figs. 4a,b the results of running the code with the initial neutron-proton 
ratio raised, for r]io = 3.0. (In these runs we have frozen the n/p ratio at the value given 
on the ordinate axis for temperatures down to Tg = 5, below which the calculation is 
allowed to proceed as normal.) Raising n/p has the effect of increasing the yields of the 
light elements, however it is important to note that in an accurate calculation the large 
yields in the high n/p regions are diluted by the smaller yields in the low n/p regions as 
well as by interactions at the interface of the regions. 

Note that increasing n/p ceases to cause significant effects for n/p^ 3, and any realistic 
calculation with multizones (e.g. Kurki-Suonio et al. 1990) has been found to have much 
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back-diffusion, reducing the magnitude of any extremes. Fig. 4b shows that in our extreme 
cases ^Be saturates at a yield of <, 10~^^, however, note that "^He is overproduced in aU 
high n/p zones by a factor of ~ 4. Thus the minimum reduction must be at least a factor 
of 4, and as mentioned before, realistic multizone models will yield even greater reductions. 
In a future paper we will investigate more thoroughly the effects of our updated network on 
inhomogeneous nucleosynthesis, however from the current preliminary exploration it seems 
unlikely that the yields can be sufficient to produce the Be and B abundances observed in 
some Pop II stars (Ryan et al 1990, Gilmore, Edvardsson & Nissen 1991, Duncan, Lambert 
& Lemke 1992). 

5. Conclusions 

We have shown that even with our extended reaction network, the standard homo- 
geneous model of primordial nucleosynthesis is unable to produce significant yields of Be 
and B. In addition, it appears that inhomogeneous nucleosynthesis is unlikely to produce 
much greater yields, even though there are uncertainties in some reactions. 

Finally, we feel it would be useful to have further data on those nuclear reaction rates 
marked S and E in table 1. In particular, ^Li(t,n)^Be, ^Li(p,7)^Be and ^Li(p,n)^Be may 
have a measurable effect on ^Be production. ^°B(a,7)^^N, ^^B(q;,7)^^N and ^^C(q;,7)^^0 
may also be critical in the production of the heavier elements. 
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Appendix 1 — Reaction Rates 
We present below the rates, in Fortran notation, used for reactions introduced or 
updated since Walker et al (1991). The notation is the usual one, used in astrophysical 
reaction rate tabulations, in which T9ab = Tg , T9Mab = Tg " and the expressions 
represent values for F = NAio'v), where Na is Avogadro's number, a the cross-section, v 
relative velocity, and the thermal average is taken over a Boltzmann distribution. 

^R + e^^Re 

F = 1.78E-9 

^Li ^ e + 24He 
F = 8.27E-1 

i^Be-^e + i^B 
F = 0.0502 

^B + e ^ 24He 
F = 9.00E-1 

12b ^e+ 12c 
F = 3.43E+1 

iiC + e^iiB 

F = 5.67E-4 

i^O + e^i^N 
F = 9.82E-3 

1^0 + e^ i^N 
F = 5.67E-3 

F = 0.0107 

F = 1.052E-4 

i^Ne + e ^ i^F 
F = 0.4146 
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^H + p^e+zz + ^H 

IF (T9 .LE. 3.) THEN 

F = 4.01E-15*T9M23*EXP(-3.380/T913) 

* (1 . + . 123*T913+1 . 09*T923+ . 938*T9) 
ELSE 

F = 1.16E-15 
END IF 

^H + e + p^zz + ^H 

IF (T9 .LE. 3.) THEN 

F = 1.36E-20*T9M76*EXP(-3.380/T913) 

* (l.-.729*T913+9.82*T923) 
ELSE 

F = 7.38E-12 
END IF 



^He + e^zz + ^H 










IF (T9 .LE. 3.) 


THEN 








F = 7.71E-12*T932*EXP(- 


.2158/T9) 




* (l.+6.48*T9+7. 


48*T9** 


2+2. 


.91*T9**3) 


ELSE 










F = 6.20E-9 










END IF 










^He + p ^ e+z/ + ^He 










IF (T9 .LE. 3.) 


THEN 









F = 8.78E-13*T9M23*EXP(-6.141/T913) 

ELSE 

F = 5.97E-15 

END IF 

'^Be + e ^ z/7 + ^Li 

IF (T9 .LE. 3.) THEN 

F = 1.34E-10/T912 * (1 . -.537*1913+3. 86*T923 

+ .0027*EXP(2.515E-3/T9)/T9) 

ELSE 

F = 6.39E-10 

END IF 

^He ^ e + ^Li 
F = 0.859 

^Li^e + ^Be 
F = 0.9846 

^Li ^ en + 2^Ile 
F = 2.9538 
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iOBe(p,7)''B 

F = 2.45E+6*T9M23*EXP(-10.39/T913) 

160(p,7)17F 

F = 1.50E+8/(T923*(l.+2.13*(l.-EXP(-0.728*T923)))) 
*EXP(-16.692/T913) 

T9A = T9/(1.+2.69*T9) 

F = 7.97E+7*T9A**(5./6.)*T9M32*EXP(-16.712/T9A**(l./3.)) 

+ 1 . 51E+8*T9M23*EXP (-16 . 712/T913) * (1 . +0 . 025*T913 

- 0.051*1923-8. 82E-3*T9) + 1 .56E+5*EXP(-6 .272/T9)/T9 

+ 1.31*T9M32*EXP(-1.961/T9) 

4He(nn, 'jfUe 

F = 4.04E-12/T9**2 * EXP(-9. 585/T9)*(l .+. 138*T9) 

^Li(n,7)8Li 

F = 3.144E+3 + 4.26E+3*T9M32*EXP(-2.576/T9) 

SLi(n,7)9Li 

F = 4.294E+4 + 6. 047E+4*T9M32*EXP(-2.866/T9) 

9Be(n,7)iOBe 

F = 1.26E+3 

i°Be(n,7)iiBe 
F = 1.32 

iiB(n,7)i2B 

F = 7.29E+2+T9M32*(2.25E+3*EXP(-0.221/T9) 

+3 . 26E+4*EXP (-4 . 514/T9) +1 . 96E+4*EXP (-10 . 804/T9) 

+3 . 90E+4*EXP (-13 . 323/T9) +5 . 86E+4*EXP (-18 . 916/T9) ) 

i60(n, 7)1^0 

F = 2.36E+l*(l.+4.45*T9)+9.66E+4*T9M32*EXP(-4.75/T9) 

1^0(11,7)180 

F = 3.11*(1.+100.*T9) 

i%(n,p)i4c 

F = 2.39E+5*(1.+.361*T912+.502*T9) 
+ 1.112E+8*EXP(-4.983/T9)/T912 

i40(n,p)i4N 

F = 2.02E+8*(1.+.658*T912+.379*T9) 
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i7F(n,p)i70 

F = 1.80E+8 

i8F(n,p)i80 

F = 1.80E+8 

F = 1.71E+11*T9M23*EXP(-10.42/T913) 

i2C(a, 7)1^0 

F = 1.04E+8*EXP(-32.120/T913-(T9/3.496)**2) 

/ (T9*(1.+.0489*T9M23))**2 

+ 1 . 76E+8*EXP (-32 . 120/T913) / (T9* ( 1 . + . 2654*T9M23) ) **2 

+ 1.25E+3*T9M32*EXP(-27.499/T9) 

+ 1.43E-2*T9**5*EXP(-15.541/T9) 

i4C(a, 7)1^0 

F = 3.375E+8*EXP(-32.513/T913)/T9**2 

+ 1 . 528E+9*T9M23*EXP (-32 . 513/T913- (T9/2 . 662) **2) 

* (l.+O. 0128*1913-0. 869*T923 

- . 0779*T9+0 . 321*T943+0 . 0732*1953) 

+ 9.29E-8*T9M32*EXP(-2.048/T9) 

+ 2 . 77E+3*EXP (-9 . 876/T9) /T9** (4 . /5 . ) 

i4N(a,7)i8F 

F = 7.78E+9*T9M23*EXP(-36.031/T913-(T9/0.881)**2) 

* (l.+O. 012*T913+1 . 45*T923+0 . 117*T9+1 . 97*T943+0 . 406*T953) 
+ 2 . 36E-10*T9M32*EXP (-2 . 798/T9) +2 . 03*T9M32*EXP (-5 . 054/T9) 
+ 1.15E+4*T9M32*EXP(-12.310/T9) 

i40(a,7)i8Ne 

F = 9.47E+8*T9M23*EXP(-39.388/T913-(T9/.717)**2) 

* (1 . + . 011*T913+1 . 974*1923+ . 146*T9+3 . 036*1943+ . 572*T953) 
+ 1.16E-1*T9M32*EXP(-11.733/T9) 

+ 3.39E+l*T9M32*EXP(-22.609/T9) 
+ 9.10E-3*T9**5*EXP(-12.159/T9) 

iiB(a,p)i4c 

F = 8.403E+15*(1.+0.022*T913+5.712*T923+0.642*T9 

+15 . 982*T943+4 . 062*T953) 

*EXP (-31 . 914/T913- (T9/0 . 3432) **2) 

+4 . 944E+6*T9** (3 . /5 . ) *EXP (-11. 26/T9) +T9M32 

* (5 . 44E-3*EXP (-2 . 868/T9) +2 . 419E+2*EXP (-5 . 147/T9) 

+4 . 899E+2*EXP (-5 . 157/T9) ) 

iiC(a,p)i4N 

T9A = T9/(l.+4.78E-2*T9+7.56E-3*T953/(l.+4.78E-2*T9)**(2./3.)) 
F = 7.15E+15*T9A56*T9M32*EXP(-31.883/T9A13) 
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i40(a,p)i^F 

F = 1.68E+13*T9M23*EXP(-39.388/T913-(T9/0.717)**2) 

* (1 . +0 . 011*1913+13 . 117*T923+0 . 971*T9+85 . 295*T943 
+ 16.061*1953) 

+ 3 . 31E+4*T9M32*EXP (-11 . 733/T9) +1 . 79E+7*T9M32 

* EXP(-22.609/T9) 

+ 9 . E+3*T9** (11 . /3 . ) *EXP (-12 . 517/T9) 

^Li(p,«7)4He 

T9A=T9/(1.+.759*T9) 

F = 1.096E+9*T9M23*EXP (-8.472/(7913)) 

-4 . 830E+8* (T9A** (5 . /6 . ) ) *T9M32*EXP (-8 . 472/ (T9A** . 333333333) ) 

+1 . 06E+10*T9M32*EXP (-30 . 442/T9) 

+1 . 56E+5*T9M23*EXP (-8 . 472/T913- (T9/1 . 696) **2) 

* ( 1 . + . 049*1913+2 . 498*1923+ . 860*T9+3 . 518*1943+3 . 080*1953) 

+1 . 55E+6*T9M32*EXP (-4 . 478/T9) 

F = 2.45E+11*T9M32*EXP(-10.39/T913) 

F = 8.57E+10*T9M23*EXP(-10.42/T913) 

iiB(p,a)24He 

F = 2.20E+12*T9M23*EXP(-12.095/T913-(T9/1.644)**2) 
* ( 1 . + . 034*1913+ . 140*1923+ . 034*T9+ . 190*1943+ . 116*1953) 
+4 . 03E+6*EXP (-1 . 734/T9) *T9M32 
+6 . 73E+9*EXP (-6 . 262/T9) *T9M32 
+3 . 88E+9*EXP (-14 . 154/T9) /T9 

i70(p,a)i4N 

F = 1.53E+7*T9M23*EXP(-16.712/T913-(T9/0.565)**2) 

* (1 . +0 . 025*1913+5 . 39*1923+0 . 940*19+13 . 5*T943+5 . 98*T953) 
+ 2.92E+6*T9*EXP(-4.247/T9) 

+ . 1* (4 . 81E+10*T9*EXP (-16 . 712/T913- (T9/0 . 04) **2) 
+ 5.05E-5*T9M32*EXP(-0.723/T9) 
+ 1.31E+1*T9M32*EXP(-1.961/T9)) 

i80(p,a)i5N 

F = 3.63E+11*T9M23*EXP(-16.729/T913-(T9/1.361)**2) 

* (1 . +0 . 025*T913+1 . 88*T923+0 . 327*T9+4 . 66*T943+2 . 06*T953) 

+ 9 . 9E-14*T9M32*EXP (-0 . 231/T9) +2 . 66E+4*T9M32*EXP (-1 . 67/T9) 
+2 . 41E+9*T9M32*EXP (-7 . 638/T9) +1 . 46E+9*EXP (-8 . 31/T9) /T9 

i8F(p,a)i50 

F = 9.52E+12*T9M23*EXP(-18.1/T913) 
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SLi(a,n)iiB 

F = T9M32*(5.505E+6*EXP(-4.410/T9) 

+4 . 596E+8*EXP (-6 . 847/T9) ) 

+1 . E+13*T9M23*EXP (-19 . 45/T913) * (2 . 02*T913 

+17 .71*1923+17 . 65*T9+3 . 57*T943) 

iOBe(a,n)i3c 

F = 3.64E+14*T9M23*EXP(-24.12/T913) 

F = 4.51E+14*T9M23*EXP(-24.33/T913) 

iOB(a,n)i3N 

F = 1.20E+13*T9M23*EXP(-27.989/T913-(T9/9.589)**2) 

iiB(a,n)i'^N 

F = 2. 468E+15*(1. +7.519*1913+1. 361*1923-14. 972*T9 

-11 . 61*1943+18 . 145*T953) *EXP (-18 . 145/T913 

-(T9/0.7207)**2) 

+1 . 459E+7*T9** (3 . /5 . ) *EXP (-11. 26/T9) +T9M32 

* ( 1 . 79*EXP (-2 . 868/T9) +1 . 678E+3*EXP (-5 . 147/T9) 

+5 . 358E+3*EXP (-5 . 157/T9) ) 

i3C(a,n)i60 

F = 6.77E+15*T9M23*EXP(-32.329/T913-(T9/1.284)**2) 

* (1 . + . 013*T913+2 . 04*T923+ . 184*T9) 
+ 3.82E+5*T9M32*EXP(-9.373/T9) 

+ 1.41E+6*T9M32*EXP(-11.873/T9) 
+ 2.00E+9*T9M32*EXP(-20.409/T9) 
+ 2.92E+9*T9M32*EXP(-29.283/T9) 

i70(n,a)i4c 

F = 3.11E+4*(1.+100.*T9)+2.12E+16*T9M23 

* EXP (-32 . 51/T913+21 . 11/T9- (2 . 33/T9) **2 . 51) /2 . 03 

i7F(n,a)i4N 

F = 7.76E+9*(l.-1.15*T912+0.365*T9)*EXP(-(T9/2.798)**2) 
+ 4.85E+10*T9M32*EXP(-15.766/T9) 

F = 6.28E+8*(1.-0.641*T912+0.108*T9) 

2H(d,7)4He 

F = 4.84E+l*T9M23*EXP(-4.258/T913) 

* (1 . + . 098*T913- . 203*T923- . 139*T9+ . 106*1943+ . 185*T953) 
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6Li(d,n)^Be 

F = 1.48E+12*T9M23*EXP(-10.135/T913) 

SLi(d,n)9Be 

F = 3.22E+11*T9M23*EXP(-10.357/T913) 

F = 4.27E+13*T9M23*EXP(-16.939/T913) 

6Li(d,p)^Li 

F = 1.48E+12*T9M23*EXP(-10.135/T913) 

^Li(d,p)8Li 

F = 8.31E+8*T9M32*EXP(-6.998/T9) 

9Be(p,d)24He 

F = 2.11E+11*T9M23*EXP(-10.359/T913 - (T9/ .520)**2) 
*(!.+. 040*T913+1 . 09*T923+ . 307*T9+3 . 21*T943+2 . 30*T953) 
+ 5.79E+8*EXP(-3.046/T9)/T9 
+ 8.50E+8*EXP(-5.800/T9)/T934 

3H(t,2n)4He 

F = 1.67E+9*T9M23*EXP(-4.872/T913) 

* (1 . + . 086*T913- . 455*T923- . 272*T9+ . 148*1943+ . 225*T953) 

^Li(t, nna)^He 

F = 8.81E+11*T9M23*EXP(-11.333/T913) 

^Li(t,n)9Be 

F = 1.46E+11*T9M23*EXP(-11.333/T913) 

9Be(t,n)iiB 

F = 3.80E+12*T9M23*EXP(-14.02/T913) 
+ 1.25E+8*T9M32*EXP(-4.43/T9) 

3He(t,d)4He 

T9A = T9 / (l.+.128*T9) 

F = 5.46E+9*T9A56*T9M32*EXP(-7.733/T9A13) 

iH(pn,p)2H 

F = 1.42E-2*T9M32*EXP(-3.720/T913) 

* (1 . + . 784*7913+ . 346*1923+ . 690*T9) 

^Be(3He, ppa)^'ile 

F = 6.11E+13*T9M23*EXP(-21.793/T913) 
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i4C(n,7)'^C 

F = 1.08E+8*3.0E-5*T9 

i3N(a,p)i60 

T9A = T9/(l.+7.76E-2*T9+2.64E-2*T953/(l.+7.76E-2*T9)**(2./3.)) 
F = 3.23E+17*T9A56*T9M32*EXP(-35.829/T9A13) 

4He(2a,7)i2c 

F = 2.79E-8*T9M32*T9M32*EXP(-4.4027/T9) 

SLi(p, na)4He 

F = 1.031E+10*T9M23*EXP(-8.429/T913) 

+6 . 79E+5*T9M32*EXP (-1 . 02/T9) 

+3 . 28E+8*T9M32*EXP (-7 . 024/T9) 

+1 . 13E+9*T9** (-0 . 433) *EXP (-3 . 982/T9) 

F = 1.08E+12*T9M23*EXP(-15.251/T913-(T9/.522)**2) 
*(!.+. 027*T913+2 . 62*T923+ . 501*T9+5 . 36*T943+2 . 60*T953) 
+1 . 19E+8*T9M32*EXP (-3 . 676/T9) 
+5 . 41E+8*EXP (-8 . 926/T9) /T912 
+4 . 72E+7*T9M32*EXP (-7 . 721/T9) 
+2 . 20E+8*T9M32*EXP (-11 . 418/T9) 

1^0(11,7)190 
F = 21.2 

20o^e + z/ + 20F 
F = 0.737 

2iF -^ e + z/ + 2iNe 
F = 0.234 



The remaining reactions are estimated as follows. 

SLi(p,7)9Be 

F = 6. 27E+8*T9M23*EXP (-8. 5/T913)*(l. +0.049*1913) 

9Li(p,a)6He 

F = 1.03E+11*T9M23*EXP(-8.533*T9M13) 

9Li(d,n)iOBe 

F = 2. 86E+11*T9M23*EXP (-10.41*7913) 

i2N(n,d)iiC 

F = 1.32E+5 
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iOBe(«,7)'^C 

F = 5.82E+14*T9M23*EX(-24.12/T913) 
+ 8.30E+7*T9M32*EX(-1.161/T9) 
+ 1.01E+8*T9M32*EX(-6.731/T9) 

F = 4.3E+15*T9M23*EXP(-6.082/T913) 
+ 6.8E+7*T9M32*EXP(-1.242/T9) 
+ 7.5E+7*T9M32*EXP(-2.832/T9) 

i2N(n,p)i2c 

F = l.OE+12 

8Li(d,p)9Li 

F = 2.58E+11*T9M23*EX(-10.34/T913) 
+ 1.24E+8*T9M32*EXP(-2.95/T9) 

9Li(p,n)9Be 

F = 1.03E+ll*T9M23*EXP(-8.533/T913) 

9Li(a,n)i2B 

F = 8.82E+15*T9M23*EXP(-19.70/T913) 

9Li(p,7)iOBe 

F = 1.03E+ll*T9M23*EXP(-8.533/T913) 
+ 3.1E+5*T9M32*EXP(-11.61/T9) 

i3N(n,7)i4N 

F = 1.32E+5 

+ 1.25E+6*T9M23*EXP(-0.16/T9) 

+ 1.32E+7*T9M23*EXP(-1.39/T9) 

iOB(a,7)i4N 

F = 5.82E+14*T9M23*EXP(-27.98/T913) 
+ 8.22E+7 *T9M32*EXP(-1.718/T9) 
+ 2.99E+10*T9M32*EXP(-1.485/T9) 

iiC(a, 7)1^0 

F = 2.05E+16*T9M23*EXP(-31.88/T913) 
+ 1.27E+8*T9M32*EXP(-0.928/T9) 
+ 1.34E+8*T9M32*EXP(-3.017/T9) 
+ 1.35E+8*T9M32*EXP(-3.365/T9) 

8B(a,7)i2N 

F = 8.67E+8*T9M23*EXP(-8.08/T913) 

* (1 . +0 . 052*T913-0 . 448*T923-0 . 165*T9+0 . 144*T943+0 . 134*T953) 
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i90(n,7)'°0 

F = 1.3E+5 

i^N(a,p)20O 

F = 3.3E+14*T9M23*EXP(-36.51/T913) 

i6C(a, 7)200 

F = 5.6E+16*T9M13*EXP(-32.82/T913) 

20F(n,7)2iF 

F = 1.3E+5 

2iF(p,a)i80 

F = 1.6E+13*T9M23*EXP(-18.10/T9) 

i7N(a,7)2iF 

F = 3.3E+14*T9M23*EXP(-36.51/T913) 
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Table 1. Reactions changed since Walker et al (1991) 



Reaction 


Notes 


Refs 


3h + e -^ 3He 


U 


T85 


8Li ^ e + 24He 


U 


T85 


i^Be^e + i^B 


A 


T85 


^B + e ^ 24He 


U 


T85 


12b -.e + 12c 


U 


T85 


iiC + e^iiB 


u 


T85 


i^O+e^i^N 


u 


T85 


1^0 +e^ i^N 


u 


T85 


i7F + e^i70 


A 


T85 


isp + e^i^o 


A 


T85 


i8Ne + e ^ 18F 


A 


T85 


iH + p^e+z/ + 2H 


A 


CF88 


^H + e + p^zz + ^H 


A 


CF88 


^He + e^zz + ^H 


A 


CF88 


3He + p^e+z/ + 4He 


A 


CF88 


■^Be + e -> z/7 + ^Li 


A 


CF88 


^He -^ e + ^Li 


A 


T85 


^Li^e + ^Be 


A 


LS78 


^Li ^ en + 2^Ile 


A 


LS78 


iOBe(p,7)iiB 


A 


W69 


160(p,7)17F 


A 


CF88 


170(p,7)18f 


A 


CF88 


4He(nn, 7)6He 


A 


CF88 


^Li(n,7)8Li 


U 


WSK89 


8Li(n,7)9Li 


A 


MF89 


9Be(n,7)iOBe 


A 


W69 


iOBe(n,7)iiBe 


A 


W69 


iiB(n,7)i2B 


US 


MF89 


i60(n, 7)1^0 


A 


W69 


i^O(n, 7)180 


A 


W69 


i^N(n,p)i4c 


U 


CF88 


i40(n,p)i4N 


A 


CF88 


i^F(n,p)i^O 


A 


W69 


i8F(n,p)i80 


A 


W69 


iiBe(p,n)iiB 


A 


W69 
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Table 1 (continued). Reactions changed since Walker et al (1991) 



Reaction 


Notes 


Refs 


i2C(a,7)^^0 


U 


CF88 


i4C(a, 7)180 


A 


CF88 


i4N(a,7)i8F 


A 


CF88 


i40(a,7)i8Ne 


A 


CF88 


iiB(a,p)i4c 


US 


CF88 


iiC(a,p)i4N 


u 


CF88 


i40(a,p)i7F 


A 


CF88 


^Li(p,a7)4He 


A 


CF88 


iOBe(p,«)^Li 


A 


W69 


iiBe(p,a)8Li 


A 


W69 


iiB(p,a)24He 


U 


CF88 


i70(p,a)i4N 


A 


CF88 


180(p,«)15n 


A 


CF88 


i8F(p,a)i50 


A 


W69 


8Li(a,n)iiB 


US 


MF89 


i°Be(a,n)i3c 


A 


W69 


iiBe(a,n)i4c 


A 


W69 


iOB(a,n)i3N 


U 


CF88 


iiB(a,n)i4N 


US 


CF88 


i3C(a,n)i60 


u 


CF88 


i70(n,a)i4c 


A 


W69 


i^F(n,a)i4N 


A 


CF88 


iSF(n,a)i5N 


A 


CF88 


2H(d,7)4He 


A 


CF88 


6Li(d,n)7Be 


A 


MF89 


8Li(d,n)9Be 


A 


MF89 


i4C(d,n)i5N 


A 


KFKM91 


6Li(d,p)^Li 


A 


MF89 


7Li(d,p)8Li 


A 


MF89 


9Be(p,d)24He 


U 


CF88 


3H(t,2n)4He 


A 


CF88 


^Li(t, nna)'^He 


A 


CF88 


^Li(t,n)9Be 


AS 


MF89 


9Be(t,n)iiB 


A 


BK89 


3He(t, dj^He 


A 


CF88 
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Table 1 (continued). Reactions changed since Walker et al (1991). 



Reaction 


Notes 


Refs 


iH(pn,p)2H 


A 


CF88 


^Be(3He,ppa)4He 


A 


CF88 


i4C(n,7)i5C 


A 


KFKM91 


i3N(a,p)i60 


U 


CF88 


4He(2a,7)i2c 


U 


CF88 


SLi(p,na)4He 


u 


BBL92 


i5N(p,a)i2C 


u 


CF88 


i80(n, 7)190 


A 


AS83 


20O ^ e + z/ + 20F 


A 


AS83 


2iF ^ e + z/ + 2iNe 


A 


EV78 


8Li(p,7)9Be 


E 




9Li(p,a)6He 


E 




9Li(d,n)iOBe 


E 




i2N(n,d)iiC 


E 




iOBe(a, 7)1^0 


E 




iiB(a,7)i5N 


E 




i2N(n,p)i2c 


E 




8Li(d,p)9Li 


E 




9Li(p,n)9Be 


E 




9Li(a,n)i2B 


E 




9Li(p,7)iOBe 


E 




i3N(n,7)i4N 


E 




i°B(«,7)14n 


E 




iiC(a,7)^^0 


E 




8B(a,7)i2N 


E 




i90(n, 7)200 


E 




17N(«,p)20O 


E 




i6C(a, 7)200 


E 




20F(n,7)2iF 


E 




2iF(p,a)i80 


E 




i^N(a,7)2iF 


E 





In the notes column, U refers to a reaction whose rate has been updated since Walker 
et al (1991); A, to a reaction which has been added; S, to a reaction with several recent 
measurements; and E, to a reaction which has been estimated. 
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The references give the most recent rate measurements and refer to: T85 — Tuh 
(1985), CF88— Caughlan and Fowler (1988), LS78— Lederer and Shirley (1978), EV78— 
Endt and Van der Leun (1978), AS83— Ajzenberg-Selove (1983), W69— Wagoner (1969), 
WSK89— Wiescher et al (1989), MF89— Malaney and Fowler (1989), KFKM91— Kawano 
et al (1991), BK89— Boyd and Kajino (1989), BBL92— Becchetti et al (1992) 
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Figure Captions 
1. The reaction network used in the code. Estimated reactions are shown with dashed 
hues. 

2a. Flow diagram for the standard model, with ryio = 3.0 (Flows to/from nuclides with 
A < 4 are neglected). 

2b. Flow diagram for the high n/p calculation, with n/p = 10 and 7710 — 3.0 (Flows 
to/from nuclides with A < 4 are neglected). 

3a. "^He mass fraction (Yp) as a function of baryon to photon ratio (ry = Uh/n^). Neutron 
lifetime is 889.6 sec. 

3b. Yields (number density relative to hydrogen) of ^H, ^He, ^Li and ^Li as functions of 
baryon to photon ratio {rj = n^/n-y). Neutron lifetime is 889.6 sec. 

3c. Yields (number density relative to hydrogen) of ^Be, ^'^B and ^^B as functions of 
baryon to photon ratio. The bands for ^Be and ^'^B are a result of the variation in 
the ^Li(t,n)^Be rate. The maximum yields for ^Be, ^'^B, ^^B within the range 2.8 < 
?7io < 4.0 are 6 x 10"^^, 2 x 10"^^, 5 x 10"^^ respectively. Within 0.01 < r]io < 100, 
maximum yields are 1 x 10~^^, 5 x 10~^^, 2 x 10~^^ respectively. 

4a. "^He mass fraction (Yp) as a function of neutron to proton ratio, for tjiq = 3.0. 

4b. Yields as a function of neutron to proton ratio, for 7710 = 3.0. 
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